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Appendix A - Excerpt from Web Page entitled, 
"Spatial Filters * Median Filter" 

found at http://www.cee.hw.ac.uk/hipr/htmynifidian.h^ 
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Median Filter 

Common Names; Median filtering, Rank filtering 

Brief Description 

The median filter is normally used to reduce noise in an image, somewhat like the mean filter . However, it often 
does a better job than the mean filter of preserving useful detail in the image. 



How It Works 

Like the mean filler, the median filter considers each pixel in the image in turn and looks at its nearby 
neighbours to decide whether or not it is representative of its surroundings. Instead of simply replacing the pixel 
value with the mean of neighbouring pixel values, it replaces it with the median of those values. The median is 
calculated by first sorting all the pixel values from the surrounding neighbourhood into numerical order and 
then replacing the pixel being considered with the middle pixel value. (If the neighbourhood under 
consideration contains an even number of pixels, the average of the two middle pixel values is used.) Figure 1 
illustrates an example calculation. 
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Neighbourhood values: 

115, 119, 120, 123, 124, 
125, 126, 127, 150 

Median value: 124 



Figure 1 Calculating the median value of a pixel neighbourhood. As can be seen the central pixel 
value of 150 is rather unrepresentative of the surrounding pixels and is replaced with the median 
value: 124. A 3 *3 square neighbourhood is used here — larger neighbourhoods will produce more 
severe smoothing. 



Guidelines for Use 

http://www.cee Jiw,ac.uk^pr/html/median,himl 6/18/2004 
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By calculating the median value of a neighbourhood rather than the mean filter, the median filter has two main 
advantages over the mean filter: 

• The median is a more robust average than the mean and so a single very unrepresentative pixel in a 
neighbourhood will not affect the median value significantly. 

♦ Since the median value must actually be the value of one of the pixels in the neighbourhood, the median 
filter does not create new unrealistic pixel values when the filter straddles an edge. For this reason the 
median filter is much better at preserving sharp edges than the mean filter. 



I shows an image dmhas been corrupted by fiaussiMLnQisa with mean 0 and s tandar d deviation (SD) 8. 

The original image is ESS for comparison. Applying a 3 *3 median filte r prod uces ElS. Note how the noise 

has been reduced at the expense of a slight degradation in image quali ty. L-ifl shows the same image with even 

more noise added (Gaussian noise with mean 0 and SD 13), and LSI is the result of 3x3 median filtering The 
median filter is sometimes not as subjectively good at dealing with large amounts of Gaussian noise as the i mean 
filter . 

Where median filtering really comes into its own is when the noise produces extreme 'outlier' pixel values, as fn 

for instance in Si which has been corrupted with 'salt and pepper 1 nni« - bits have been flipped with H 

probability 1%. Median filtering this with a 3 *3 neighbourhood produces ES, in which the noise has been $ 
entirely eliminated with almost no degradation to the underiying image. Compare this with the similar test on =- 
the mean filter. £ 

W 

Consider another example wherein the onggial unage 120 has been corrupted with higher levels (i.e. p=5% that HI 

a bit is flipped) of salt and pspprr noHf* r™. After smoothing with a 3 x3 filter, most of the noise has been O 

eliminated E2E1. If we smooth the noisy image with a larger median filter, eg 7*7, all the noisy pixels ^ 

disappear, as shown in QH Note that the image is beginning to look a bit 'blotchy 1 , as greylevel regions are 
mapped together. Alternatively, we can pass a 3x3 median filter over the image three times in order to remove 

all the noise with less loss of detail uSj. 

In general, the median filter allows a great deal of high spatial frequency detail to pass while remaining very 
effective at removing noise on images where less than half of the pixels in a smoothing neighbourhood have 
been effected (As a consequence of this, median filtering can be less effective at removing noise from images 

corrupted With flaussifln noi<ifl ) e 

One of the major problems with the median filter is that it is relatively expensive and complex to compute. To 
find the median it is necessary to sort all the values in the neighbourhood into numerical order and this is 
relatively slow, even with fast sorting algorithms such as quicksort. The basic algorithm can however be 
enhanced somewhat for speed. A common technique is to notice that when the neighbourhood window is slid 

htto://ww.cee.hw.ac.uk^ 6/1 g/2004 
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across the image, many of the pixels in the window are the same from one step to the next, and the relative 
ordering of these with each other will obviously not have changed. Clever algorithms make use of this to 
improve performance. 

Exercises 



mm 



1. Using the image JwflB , explore the effect of median filtering with different neighbourhood sizes. 

2. Compare the relative speed of mean and median filters vising the same sized neighbourhood and image. 
How does the performance of each scale with size of image and size of neighbourhood? 

3. Unlike the mean filter , the median filter is non-linear. This means that for two images A(x) and B(x) : 
median[A(x) + B(s]) ^ fr^dian[A{x)] + median[B{x)] 

Illustrate this to yourself by performing smoothing and pixel addition (in the order indicated on each side 
ofthe aboveejraation!) to a set of test images. Carry out this experiment on some simple images, e.g. 

landL— I 
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Local Information m 

O 

General advice about the local HIPR installation is available hejs § 
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Hypermedia Image Processing Reference 
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